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ABSTRACT 

We show that using either the method of Page & Carrera or the well-known 1 /V a method to construct the 
binned luminosity function (LF) of a flux limited sample of Active Galactic Nuclei (AGN) can produce an 
artificial flattening (or steepening in the case of negative evolution) of the binned LF for bins intersected by the 
flux cutoff of the sample. This effect is more pronounced for samples with steep and strongly evolving parent 
LFs but is still present even for non-evolving LFs. As a result of this distortion of the true LF, fitting a model 
LF to binned data may lead to errors in the estimation of the parameters and may even prompt the erroneous 
use of broken power law functions. We compute the expected positions of apparent breaks in the binned LF. 
We show that these spurious breaks in the binned LFs can be avoided if the binning is done in the flux-redshift 
plane instead of the typically used luminosity-redshift plane. Binning in the flux-redshift plane can be used 
in conjunction with the binning in the luminosity-redshift plane to test for real breaks in the binned LFs and 
to identify the features that are the result of binning biases. We illustrate this effect for most typical forms of 
luminosity dependence and redshift evolution and show how the proposed method helps address this problem. 
We also apply this method to the MOJAVE AGN sample and show that it eliminates an apparent break in the 
binned LF. 

Subject headings: galaxies : luminosity function, mass function — galaxies : evolution — galaxies : active — 
quasars : general — 



1. INTRODUCTION 

The luminosity function (LF) of active galactic nuclei 
(AGN) and its redshift dependence is highly useful for study- 
ing the cosmological evolution of AGNs and for statistical 
tests of AGN unification models. Typically the binned LF 
is obtained using the standard 1 /V max method (Schmidt 1968; 
iFelten II 19761) or its gen eralization for combined samples 1 /V a 
(lAvni & Bahcal II 119801) . The evolution of the LF is estimated 
using the (V/V max ) method (ISchmidtl fr968). Sometimes a 
model L F is derived by fitting a model LF to the binned LF 
(see, e.g.. lPadovani et aLll2007l) . However, the binned LF ob- 
tained using the 1 /V a method suffers from se veral biases (see 
lLa Franca & Cristianil 1 1 9971: IWisotzkilll998l and references 
therein) and for this reason many prefer fitting the model LF to 
the unbinned data with, e.g., maximum likelihood techniques. 
In this case the binned LF, however, is often used before per- 
forming the fit to observe the overall behavior of the luminos- 
ity function, e.g., to decide its evolution type and to decide 
whether or not a double power law LF should be fitted to the 
data instead of a single power law. It can also be used after 
th e fit to visualize t he goo dness of fit. 

Page & Carrera (2000) described a different method of 
computing the binned LF that, according to the authors, gives 
better estimates of the LF for lower luminosity bins that are 
close to the flux limit of the survey using an improved (over 
1 /V a ) estimator. However, their method does not eliminate all 
the biases associated with undersam pled bins that are n ear the 
flux cutoff o f the survey (see, e.g ., lM~ivaii et al.ll2001l) . In a 
recent paper (ICara & Lister 2008) we derived the radio lumi- 
nosity function of the M OJAVE survey of core-selected AGN 
dLister & H oman 2005). While a single power law provided 



a good fit to the u nbinned data, the binn ed LF, computed us- 
ing the method of lPage & Carrera! (|20Q0), was flatter at lower 
luminosities, and sugg ested a double power law LF (see, e.g., 
l Arshakian et al.ll2.Q06l) . After an analysis of the problem (see 
ICara&Listei| [2008) we concluded that the flattening of the 
binned LF was due to the fact that the values of the LF com- 
puted over small chunks of undersampled (because of the flux 
cutoff) bins are not good approximations of the value of the 
LF at the centers of the bins, particularly when the LF is a 
steep and strongly evolving function across the bin. 

There are many other examples in the literature where 
binned luminosity functions have been e xtended well below 
the flu x cutoff of the survey. For example, Dunlo p~& Peacock! 
(1990, Fig. 10) used the binned LF to compare the predic- 
tion of their pure luminosity evolution (PLE) models with the 
data. The plot displays a clear flattening of the LF for bins 
below the flux cutoff. Because of the multiple flux limits in 
their data, this flattening happens gradually and, therefore, the 
binned LFs in that plot are better described as being "bent" 
than "broken". 

A somewhat dif ferent example of using bins below the flux 
cutoff comes from lUrrv & Padovanild 1995b . In Figure 14 the 
authors present the local LFs of the FR II galaxies, steep spec- 
trum radio quasars (SSRQ) and flat spectrum radio quasars 
(FSRQ). In this case, because the luminosities have been de- 
evolved, the LFs have been constructed essentially using one 
large redshift bin with a modified (de-evolved) flux limit. 
Therefor e, almost all the lumino sity bins presented in Fig- 
ure 14 of lUrrv & Padovanil (1 19951) lie below the survey limit 
(however, the precise number of bins below the limit is dif- 
ficult to determine because the redshift bin width is not pre- 
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sented). At very low luminosities, particularly at luminosities 
where the survey limit starts cutting into the number of objects 
in a bin, the number densities will start decreasing because the 
incomplete bins, due to the shape of the flux cutoff, will pref- 
erentially sample a different part of the LF than would be the 
case for a complete bin. We will show in Section |2] that the 
position of the dominant break will depend on the minimum 
redshift of the bin in the case of a non-evolving luminosity 
function. 

A final example is that of the binned luminosity function of 
the quarter-Jansky sample from lWall et al.l d2005l) . The flat- 
tening of the LF of this sample at low luminosities produces 
the illusion of the luminosity evolution, which was recognized 
by the authors as such. The authors attributed the flattening 
of the LF at lower luminosities to an intrinsic spread in radio 
spectral indices. 

Since this flattening (or steepening in the case of nega- 
tive evolution) is caused by undersampling of the bins cut 
by the survey limit, one solution to eliminate these spuri- 
ous breaks is to stop binning as soon as the lowest lumi- 
nosity bin reaches the survey limit. This, however, can lead 
to the exclusion of data points lying between this last com - 
plete bin and the su rvey limit. La Franca & Cris tianil (119971) : 
Mivai i et al.l d2001l) proposed to correct for biases in the un- 
dersampled bins by prorating them for the missing parts us- 
ing a model LF determined by independent means. This ap- 
proach, howev er, has the disadvantag e of depending on an as- 
sumed model. Efsta thiou et al.l d 1 988b have developed a non- 
parametric "step-wise maximum likelihood" method that al- 
lows a binned LF to be fit directly to the unbinned data. While 
these methods their own advantages and disadvantages, we 
have chosen to focus in this paper on the method of Page & 
Carrera. We note that any method that computes the values of 
the LF over incomplete bins must, in one way or another, use 
some sort of extrapolation to the parts of the bins lying below 
the flux limit. 

In this paper we show that the biases associated with un- 
dersampled bins can produce artificial breaks in the derived 
luminosity functions, particularly in the case of strong evolu- 
tion. As evolution of the LF is typically judged by the shifts 
of features in the LF, these artificial breaks, when interpreted 
as real, can produce the illusion of luminosity evolution. We 
show that the flattening (or breaking) of the binned LF for 
bins close to the flux cutoff of the survey can be minimized 
if the binning is performed in the flux-redshift plane instead 
of the luminosity-redshift plane as it is usually done. The im- 
provement is most apparent at lower redshifts where, due to 
the steepness of the luminosity cutoff in the logL-z space, 
many bins are undersampled, and therefore give strongly bi- 
ased estimates of the luminosity function. While our analysis 
of the distortion of t he shape of the binned LFs is performed 
using the method of Pag e & Carrera! d2000t) for constructing 
binned luminosity functions, this analysis equally applies to 
the l/y a method as well. The binned LFs constructed using 
the 1/V a method will e xhibit larger dev i ations from the true 
LF than the method of iPage & Carrera! d2000t) at extremely 
low luminosities (L < Lbieak,3, see equation |9}. However, in 
the case of non-evolving LFs, the 1 /V a method, unlike the 
method of Page & Carrera, will not distort the shape of the LF 
for Lbreak,3 < L < Lb rea k,i because it assigns a variable weight 
to each object in the bin and can correctly handle the flux 
limit down to lower luminosities (Lbreak.3)- Regardless of the 
method used, any determination of the LF below Lbieak,3 is 
pure extrapolation and should be avoided. 



The outline of the paper is as follows: In § [2] we de- 
scribe both binning methods and discuss the main advantages 
of the binning in the flux-redshift plane over the traditional 
luminosity-redshift plane binning method. In § [3] we apply 
these two binning methods to several model luminosity func- 
tions and show that the LF binned in the flux-redshift plane 
is a better approximation to the model LF. We construct the 
binned LF of the MOJAVE sample using both methods in §[4] 
We summarize our findings and present our conclusions in 
§0 

Throughout this paper we assume a cosmology with fi m = 
0.3, 0\ = 0.7, Q, r = and H = TOkms-'Mpc" 1 . All lumi- 
nosities are quoted as monochromatic luminosities at specific 
frequency v. We also adopt the following convention for the 
spectral index, a lat j: S u oc v ami and assume that a, a( j = for 
all model LFs and plots in this paper. 

2. METHOD 

The differential luminosity function of a sample of astro- 
physical objects is defined as the number of objects per unit 
comoving volume per unit luminosity interval: 



d 2 N 



(1) 



dVdL 

but, because of the typically large span of the luminosities, it 
is often defined in terms of logL: 



(j)(\ogL,z) ■■ 



d 2 N 



(2) 



dVd log L 

and the relationship between the two forms is given by: 

(j) (logL,z) = In 10 L4>(L,z). (3) 

One of the advantages of working with the LF defined in terms 
of logL is that for the most common forms of the LF used in 
model fitting, i.e., the simple power law 
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the plot of \ogcj) vs. logL will show a simple linear depen- 
dence: 

log </> (logL, z) = (a+ l)logL+log(lnlO)-(a + l)logL* (5) 

with a slope (a+ 1). In equations (0| and © L* is an arbitrary 
positive constant. 

The usual way to construct the binned LF of a sample of 
objects is to divide the logL-z plane into luminosity-redshift 
bins and to estimate the space density in a bin using, for ex- 
ample, the 1 /V max method or other methods. In thi s paper we 
use the improved method of Page & Carrera! d2000l) . that is, 
we assume that the LF at the center of a bin with a luminosity 
interval logL m j n and logL max and a redshift interval z m i n and 
Zmax can be approximated as: 



</>(logL c ,Z c ): 



N 



fZmax (-logLmax dV J J 1 T ' 



(6) 



where N is the number of sources detected within the bin lim- 
its and L m j n (z) is the minimum luminosity within the bin at 
which we can still detect an object. In Figure Q] we show a 
luminosity-redshift bin (hatched bin "A") positioned so that 
it just touches the flux cutoff of the survey at z = 0.5. Any 
bin of the same luminosity width A logL sampled in the 
logL - z plane and centered below the center of the bin "A" 
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will have some parts of it lying below the survey's flux limit. 
If the luminosity function of the parent sample is steep and 
strongly evolving across the bin, then equation © will pro- 
duce stronger biases for undersampled bins than for bins com- 
pletely above the flux cutoff. We note that the shape of the 
binned luminosity function will not be affected by biases as 
long as these biases are constant for all bins in a given red- 
shift interval. However, this is not the case for the undersam- 
pled bins: given a negative power law index for the LF in the 
parent sample (a < 0), in the case of positive evolution, the 
change in bias will tend to flatten the LF, while in the case of 
negative evolution, it will tend to steepen the LF at luminosi- 
ties near and below the luminosity limit. The position of this 
break, using L = 4t:SDI(z)(1 + z )~ il+a "' 1 \ will be given by 



log Lb, 



reak. 1 



log 



47r5 m i n £) L (Zmax) 
(l+Zmax) 1+ "'» d 



+ AlogL/2, (7) 



where S m i n is the flux limit of the survey and £>l(z) is the lumi- 
nosity distance. This break luminosity corresponds to the lu- 
minosity coordinate of the center of the bin "A" shown in Fig- 
ure!]] The break (i.e., the difference between the slopes before 
and after the break luminosity) will be more pronounced for 
more strongly evolving LFs. Another change in shape of the 
LF is possible when the flux cutoff line passes through the 
upper-right (high redshift and high luminosity) corner of the 
bin: 

logLbreak.2 = logL bre ak,l ~ AlfJgL. (8) 

Finally, the last break will occur around 



log^breakJ =log 



47rS m ; n Z)^(z m i n ) 



+ AlogL/2, (9) 



(1+Zmm) 1+Q '" d 

which is the coordinate of the center of the bin whose lower- 
left (low redshift and low luminosity) corner is just touching 
the survey limit. For non-evolving LFs, this is the dominant 
break, since at luminosities lower than Lbreak,3 we will be sam- 
pling only the lower density part of the bin. The coordinates 
of the last two breaks are approximate because their effec- 
tive position will depend on the bin dimensions, steepness of 
the survey limit within the bin limits and steepness of the LF 
itself. In addition, if a rac j is not constant, the positions of the 
above breaks will be blurred since the breaks will occur at dif- 
ferent luminosities for sources with different spectral indexes. 
Alternatively, one can think that each source has it s own sur- 
vey li mit (similar to the "single source survey" in iWall etal.1 
2005) and the observed break will be the result of overlapping 
of all these individual breaks. 

In the case of samples with multiple limits, the breaks given 
by equations (0, (© and (0 should be computed for each 
flux limit. That is, the total number of breaks in the binned 
(in logL-z plane) LF potentially will be equal to three times 
the number of flux limits. The strengths of each break will 
depend on the sky area covered for each flux limit (i.e., by the 
contribution of each subsample). The first break, however, 
will occur at L = Lbreak.i computed for the largest flux limit in 
the sample. 

The undersampled bins distort the shape of the LF and in 
principle should not be used in the analysis. Doing this, how- 
ever, would discard a significant amount of information since, 
in samples with a < 0, sources tend to accumulate closer to 
the luminosity limit of the survey. The loss of information, 
however, can be avoided if the sampling is done in the flux- 
redshift plane while keeping all bins above the flux cutoff and 
thus avoiding the distortion of the LF shape. Sampling in 
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Fig. 1. — A typical bin ("A") in logL-z space is shown together 
with a bin ("B") that is a rectangle in logS-z space. Both rectangles 
lie just above the flux limit of the survey (solid line L m i n (S m i„,z)). 
The bins' height is AlogL = AlogS = 0.5 and they span a red- 
shift range 0.04 < z < 0.5. The dash-dotted line shows the upper 
limit of the bin in the logS-z plane (line Lmin(5'min,z) + AlogL). 
The bin "A" is centered on logLb re ak,i ~ 27.23 and spans a range 
of luminosities 26.98 < logL < 27.48. The bin corresponding to 
logLbrcak.2 ~ 26.73 that would span a range of luminosities 26.48 < 
logL < 26.98 (i.e., it would be adjacent to bin "A") and the bin cor- 
responding to logLbreak,3 ~ 24.98 that would span a range of lumi- 
nosities 24.73 < logL < 25.23 are omitted here for clarity, 
the log S — z plane was used previously by Koch anekl (1 1996b 
to minimize the effects of cosmology and variations in the 
spectral index. Here we will convert the binned LF back to 
logL-z coordinates to compare the shapes of the LFs ob- 
tained through both methods. The method of constructing the 
binned LF in the logS-z space is essentially similar to the 
method described above. Noting that for constant spectral in- 
dexes dlogS = if logL and that the limits of integration in the 
equation ^ are fixed in the logS-z space, the double integral 
can be easily computed and the value of the binned LF can be 
estimated as 



0(logS c ,Zc) ' 



N 



(Vfonax)-V(Zmm))AlogS' 



(10) 



where z m ; n and z max are the redshift limits of the bin, and V(z) 
is comoving volume. 

We have mentioned above that the distortion of the LF is 
strongest at undersampled bins, but even the complete bins 
are not fully immune from distortions. However, for complete 
bins, we expect the bias, i.e., the difference between the true 
value of the LF at the center of the bin and the value computed 
using either equation © or ( [Tol l, to change little, at least for 
well behaved parent LFs such as the single power law LF from 
equation (01 or a smoothed two power law LF of the form 



<j>(L,z) = 



4>o/L hrt 



y^-break / \^break 



(11) 



at least for pure density evolution (PDE). In the case of pure 
luminosity evolution (PLE) the two power law LF may exhibit 
more complex biases and possible larger distortions, but these 
should be confined mainly to the region around the break lu- 
minosity L break . 

From Figure Q] we see that the complete bins in the two 
methods, when centered around the same luminosity, sample 
different regions of the logL-z plane and even if the binned 
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LF may be distorted, we would expect these distortions to be 
different. This can be used to test whether a break seen in a 
binned LF is real: a real break should appear with high prob- 
ability in the binned LFs regardless of the sampling method, 
while a spurious break that is due to binning biases will de- 
pend on the sampling method. 

Another immediate observation from Figure Q] is that the 
centers of the complete bins closest to the luminosity limit 
are lower for bins sampled in the \ogS-z plane. It can be 
shown that for the lowest bins in the two sampling methods 
(bins "A" and "B" in Figure[T]l 



log^=log^-(a rad+ l)logfi±^f: 

V i+z, 



-Az/2 \ 

nax / 



+2 log | 



■Az/2) 



Dh (Zmax) 



(12) 



so that logLj 3 is smaller than logL^ as long as D\iz){\ + 
^-(i+Qnid) j s an increasing function of redshift. Thus, if we limit 
ourselves to complete bins only, the sampling in the logS-z 
plane allows the LF to be determined down to lower luminosi- 
ties than when sampling in the logL-z plane. 

Finally, for steep LFs, the density of the sources will be 
greater along the flux limit of the survey so that the lower lu- 
minosity bins sampled in the log S— z plane will contain more 
sources than the bins sampled in the logL-z plane. This al- 
lows one to reduce the luminosity width of the bins and obtain 
a finer sampling of the LF near the flux limit. 

One disadvantage of the sampling in the logS-z plane is 
that the bins, when converted to the logL-z plane, span a 
larger range of luminosities when compared to direct sam- 
pling in the logL-z plane. This can result in blurring of real 
features (if present) in the binned LF. This blurring can be 
somewhat minimized by reducing the redshift width of the 
bins at the expense of higher Poisson errors in the high lumi- 
nosity bins. 

3. COMPARISON OF BINNED AND MODEL LUMINOSITY 
FUNCTIONS 

In this section we construct the binned LF for several com- 
monly used forms of model luminosity functions and their 
evolution. We compute the binned LF in the logL-z plane 
from equation §6^ with N computed as: 

AWi= / / 0(logL,z)— dzdlogL. (13) 

Similarly, for bins sampled in the \ogS-z plane we use 

N mode i = / / 0(logL(5,z),z)— dzdlogS. (14) 

-'Zmin -nosrSmin 



dz 



Alternatively, one could use Monte Carlo simulations to pro- 
duce the desired samples and construct the binned LFs from 
these simulated samples. However, these samples would need 
to be very large to minimize the low number statistics at the 
extreme luminosity ends of the binned LFs. In this section we 
will assume a survey flux limit S m j n = 1.5 Jy. The parent pop- 
ulation has redshift range z\ = 0.05 and zi = 4, and luminosity 
range from U = 5 x 10 24 WHz"' to U = 10 30 WHz" 1 . These 
values are based on lCara & Lister! d2008l) . We will normalize 
the LFs so that N = 1. 

The biases in the binned LFs will be different for the two 
sampling methods, but since we want to compare the distor- 
tions of the LFs, a constant shift of the LF is not important for 
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Fig. 2. — Differential (dotted line) and binned luminosity func- 
tions in logL-z plane using the method of Page & Carrera (dashed 
line) and in logS-z plane (continuous line) for a non-evolving 
single power law LF with a = -2.5. The redshift intervals are 
0.05 < z < 0.4, 0.4 < z < 1, 1 < z < 2.2 and 2.2 < z < 4. The LFs 
in higher redshift intervals are shifted upwards along the density axis 
by 1 dex compared to the LF in the previous redshift interval. The 
thin vertical lines show the position of the last complete bin (Lbieak.i) 
in the logL-z plane and the thick vertical lines show the position of 
the break Lb re ak,3- The solid lines are truncated at luminosities given 
by equation dl2t with logL c = Lbreak.i, below which the LFs binned 
in logS-z plane would be extrapolated. The binned LFs were com- 
puted using AlogL = 0.5. 

our purposes. Therefore, in order to make it easier to compare 
the model LF with the binned LFs we shift the binned LFs by 
an amount equal to the bias of the first (i.e., lowest luminosity) 
complete bin, computed as (4>} bin - center 

We begin by showing the effects of the binning in logL-z 
plane on a moderately steep (a = -1.5) single power law LF 
with no evolution. The binning is performed in four redshift 
intervals with AlogL = 0.5. The results are presented in Fig- 
ure |2] with the LFs in higher redshift bins shifted upward by 
1 dex compared to the LF in the previous bin. The position 
of the beginning of the break corresponding to Lbreak.i from 
equation dT) is marked by a thin vertical line for every LF and 
the thick vertical lines show the positions of Lb rea k.3- The ends 
of the solid lines correspond to logL^ from equation ( fl2b with 
logL^ = Lbreak.i for each redshift bin. Although LFs should 
never be calculated below Lbreak,3, we plot them here only to 
illustrate the significant deviations that result. Several obser- 
vations can be drawn from Figure 12 1) even for a moderately 
steep non-evolving LF, binning below the survey limit in the 
logL-z plane produces a broken LF; 2) the effective posi- 
tion of the break will depend on the redshift limits of the bins, 
and for a non-evolving LF the break is closer to Lbreak.3 from 
equation (O; 3) the break shifts along the luminosity axis, 
thus mimicking luminosity evolution; 4) the LF binned in the 
logS-z plane preserves the shape of the model differential 
LF. While the difference between the two binned LFs is neg- 
ligible over the common luminosity range, it becomes more 
pronounced for steeper and stronger evolving functions. To il- 
lustrate the effect of steepening of the LF, in Figure [3] we plot 
the binned LFs in the first redshift interval (0.05 < z < 0.4) 
for three values of the model LF slope: a = -1 .5, a = -2, and 
a = -2.5. As expected, the break is more evident for steeper 
LFs. We note here that for non-evolving LFs, the 1 /V a method 
will correctly reproduce the shape of the LF down to Lbreak.3 
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Fig. 3. — Differential (dotted line) and binned luminosity functions 
in logL-z plane (dashed line) and in logS-z plane (continuous line) 
for a single power law LF for three values of the slope a = -1.5, 
a = -2 (shifted upwards by 1 dex compared to the a = -1.5 LF) and 
q = -2.5 (shifted upwards by 3 dex compared to the a = -1.5 LF). 
The redshift interval is 0.05 < z < 0.4. The thin vertical lines show 
the position of the last complete bin (Lb rea k.i) in the logL-z plane 
and the thick vertical lines show the position of the break Lb re ak,3- 
The binned LFs were computed using AlogL = 0.5. 

because its variable weight can correctly handle the flux cut- 
off for bins above Lbreak.3, while the method of Page & Carrera 
assumes a constant value of the LF within a bin. 

We now explore the effects of binning on LFs containing 
pure density evolution. There are several analytic forms of 
evolution often used in modeling of the luminosity functions: 
(1 +z) , exp[fcr(z)] (where r(z) is the look-back time in units 

of Hubble time) and exp — i (^ a ) 2 . The difference between 

the first two evolution forms is minor for our purposes and we 
choose the second form of evolution with k = 5. This value of 
the evolution parameter i s consistent with the value reported 
by iPadovani etail d2007l) for the DXRBS flat spectrum ra- 
dio quasar sample when the latter is converted from PLE to 
PDE. The binned and differential LFs with exponential den- 
sity evolution and slope a = -2.5 are presented in Figure [4] 
Compared to the non-evolving LF, the break in the LF binned 
in the logL-z plane is now much more pronounced, while 
the LF binned in the log S- z plane does not show any break. 
Intuitively, this is because the mean redshifts in undersam- 
pled bins deviate from the mean redshifts in complete bins by 
larger values when evolution is present. Again, the shift of 
this break along the luminosity axis produces the illusion of 
luminosity evolution. 

The Gaussian form of evolution is different from the first 
two forms because it produces strong negative evolution for 
Z > zo- For this form of evolution we u se the parameter val - 
ues for model B in cosmology III from 1 Willott et al.l (| 1998): 
Zo = 2.34 and a = 0.631 with luminosity slope a = -2.5. The 
results are presented in Figure [5] Compared to the previous 
evolution model in Figure [4] the apparent breaks are much 
stronger in the lower redshift bins, however the break almost 
disappears in the highest redshift interval. This is due to nega- 
tive evolution for z > 2.34, which actually increases the slope 
of the LF in the lower luminosity range. For this particular 
selection of parameters, the steepening of the LF due to the 
negative evolution is balanced somewhat by the artificial flat- 
tening created by the binning method. To better illustrate this 



Fig. 4. — Same as in Figure [2] except the LF has pure density 
evolution (PDE) of the form exp [kr(z)] with k = 5. The luminosity 
slope is a = -2.5. 
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Fig. 5. — Same as in Figure [2] except the LF is evolving with PDE 
of the form exp [-1 (^) 2 ] with zo = 2.34 and a = 0.631. The lu- 
minosity slope is a = -2.5 for the lower four LFs and a = -1.5 for 
the upper LF. The LFs in higher redshift intervals a shifted upwards 
along the density axis by 1 dex compared to the LF in the previous 
redshift interval and the LF in the 2.2 < z < 4 with a = -1.5 was 
shifted upwards one more dex compared to the LF with a = -2.5 
slope in the same redshift interval. The thin vertical lines show the 
position of the last complete bin (Lb rea k,i) in the logL-z plane and 
the thick vertical lines show the position of the break Lbreak,3- 

point we replotted the LF in the last redshift interval using a 
luminosity distribution slope a = -1 .5. In this case the flatten- 
ing due to steep luminosity dependence cannot balance com- 
pletely the steepening due to negative evolution. 

If the artificial offsets are removed from the LF curves in 
Figure [5] one can see that unlike the case o f the bias due to 
spectral spread reported bv lWall et al.l (|2005), the apparent lu- 
minosity evolution due to spurious breaks in the binned LFs 
is created by the breaks at different space densities. This is 
in good agreement with decreas e in the density at which the 
breaks occur seen in Figure 3 of IWall et al.l J2005l) . Also, the 
steepening below the breaks of the LF in the case of nega- 
tive evolution is in good agreement with the steepening of 
the binned LF in the highest redshift interval in Figure 3 of 
IWall et all d2005f) . The negative evolution in that redshift in- 
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Fig. 6. — Luminosity functions for a two power law luminosity dis- 
tribution with ai = -1.5 and aj = -3 and a sharp break. We assume 
exponential PDE of the form exp[A:r(z)] with k = 5. Notations are 
the same as in Figure(2] In this plot only the LF in the 2.2 < z < 4 
redshift interval was shifted upwards by 1 dex. 

terval is revealed by the shift of this function to the left (to- 
wards lower luminosities). We point out that the break in the 
LF due to undersampled bins (Lbreak.i) will occur at larger lu- 
minosities than the ones corr esponding to the vertical bars in 
Figure 3of lWall etaHd2005l) . since their sample has multiple 
flux limits. As discussed in Section[2] the flattening will begin 
at Lbreak.i, computed for the largest flux limit of the survey. 

Finally, we examine how real features in LFs are blurred by 
different sampling methods. For this purpose we will use a 
two power law model LF with a sharp break: 

4>(L z) = 4>o { (^/^ break ) '' Ll < L < ^ break ' (15) 

1 (i/ibreak) 2 , ^break < L < L2, 

with a\ = -1 .5, ct2 = -3, and logLbreak = 28 with exp [5r(z)] as 
the parametrized density evolution. The results are presented 
in Figure [6] Again, the LFs binned in the logL-z plane di- 
verge from the differential LF below the luminosity limit and 
yield a flatter low luminosity slope. Unlike the single power 
law cases discussed previously, in this case the real break in 
the LF eliminates the appearance of luminosity evolution. We 
can see that with binning in the logS-z plane it is possible 
to reach lower luminosities than logL-z binning without pro- 
ducing artificial breaks. However, this comes at the expense 
of blurring the real breaks: we can see that the (real) break 
in the LF binned in the logS-z plane is wider than the (real) 
break of the LF sampled in the logL-z plane. This blurring 
can be minimized somewhat by reducing the redshift width of 
the bins at the expense of smaller number statistics in the high 
luminosity bins. 

4. THE BINNED LF OF THE MOIAVE SAMPLE 

As an illustration using actual data, we construct the 
binn ed LF of the MOJAVE sample for the FR Il-only sam- 
ple dLister & Homanl 120051) using binning in the logS-z 
plane and binn i ng in logL-z plane using the method of 
Page & Carrera (2000). The LFs constructed in t he log L-z 
plane were presented in Figure 4 of Car a & Lister! (120081) . Be- 
cause the distortions of the LF are predominant in the red- 
shift intervals of strong evolution, and because the only red- 
shift interval in which there are enough LF points to make 
any conclusions about the distortions induced by binning is 
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Fig . 7 . — MOIAVE LF of FR Il-only sources for the redshift interval 
0.04 < z < 0.7 and AlogL = 0.5 (upper panel) and AlogL = 0.25 
(lower panel). The continuous line shows the differential model LF, 
while the dashed line shows the binned model LF. The filled circles 
show the observed LF binned in the logL-z plane and open squares 
- the observed LF binned in the logS-z plane. The gray arrows 
indicate the position of the lowest luminosity complete bin in the 
logL-z plane corresponding to Lbreak.i • 

0.04 < z < 0.7, we choose this redshift interval for our il- 
lustration. In Figure [7] we present the results of construct- 
ing the binned LF using the two methods of sampling with 
AlogL = 0.5 (upper panel) and AlogL = 0.25 (lower panel). 
The continuous line shows the differential model LF derived 
using a maximum likelihood method (f or a description of th e 
procedures and model parameters, see Car a & Lister! 2008). 
The dashed lin e shows the binned mo del LF computed using 
the method of iPage & Carreral (120001) in the logL-z plane. 
In order to facilitate the comparison of the distortions of the 
shape of the LFs we have compensated for the constant bias 
as described in Section[3] that is, we shifted the binned LFs by 
an amount equal to the bias of the lowest luminosity complete 
bin computed as (4>) hin - 01 bin center- The observed LF binned 
in the logL-z plane and represented by filled circles agrees 
w ell with the binned mod el LF constructed using the method 
of Pag e & Carreral (|2000), that is, the model provides a good 
fit to the data. However, both these luminosity functions flat- 
ten below Lbreak.ii producing an apparent break. On the other 
hand, the observed binned LF constructed in the log S-z plane 
shows no sign of flattening for L < Lb re ak,i- It also confirms 
our conclusion that the apparent break in the binned MOJAVE 
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LF is due to strong evolution across undersampled bins. The 
plot in the lower panel was created using a smaller luminos- 
ity bin size (AlogL = 0.25). One can see that the LF binned 
in the logS-z is more tolerant to smaller bin sizes near the 
flux limit. This is because for luminosity functions with steep 
(negative) luminosity slopes, the sources will tend to accumu- 
late near the survey limit. 

5. SUMMARY AND CONCLUSIONS 

The binned LFs constructed in the traditional logL-z plane 
show significant distortions in bins that are undersampled due 
to luminosity limits of the survey. In the presence of strong 
positive evolution across these bins, the derived luminosity 
function will flatten for luminosities below the break lumi- 
nosities given by equations ||7), (© and ||9), producing artifi- 
cial breaks in the binned LFs. In the case of negative evolution 
the effect is opposite: the luminosity functions steepen. For 
single power law LFs such breaks can create the illusion of lu- 
minosity evolution. The space density at which these spurious 
breaks occur varies with the redshift bin. 

In order to avoid the possibility of producing spurious 
breaks in binned LFs, the binning in the logL-z plane should 
not extend below luminosities given by Lbreak.i from equa- 
tion (0 and binning below Lbreak,3 should be avoided at all 
costs. However, this may mean that some sources are ex- 
cluded from the analysis. Alternatively, we propose that the 
binned LFs be constructed in the logS-z plane using com- 
plete bins only. This method has the following advantages: 

1. Sampling in the logS-z plane can be done using com- 



plete bins without the loss of information (i.e., exclu- 
sion of sources with L < Lbreak.i) that would happen if 
only complete bins were used in the logL-z plane. 

2. Using complete bins eliminates the appearance of spu- 
rious breaks or steepening of the LFs. 

3. In the same redshift interval the LFs binned in the 
log S— z plane can sample the parent LF down to lower 
luminosities than in the case of logL-z binning as long 
as D\{z)(l +z)" (1+arad) is an increasing function of red- 
shift (see equation [T2l. 



4. The previous advantage comes at the expense of blur- 
ring the real features in the binned LF. 

5. Binning in the \ogS-z plane should be used in conjunc- 
tion with binning in the logL-z plane (for L > Lbreak,i) 
to test for real breaks in the binned LFs and to identify 
the features that are the result of binning biases, since 
the biases associated with binning in the two planes are 
expected to be different. This is especially important 
because the binning in logS-z may blur some weak 
real features in the binned LF. 



We thank the anonymous referee for critical comments that 
have improved this paper. 
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